Edge Excitations and Non-Abelian Statistics in the Moore-Read State: A Numerical 
Study in the Presence of Coulomb Interaction and Edge Confinement 
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We study the ground state and low-energy excitations of fractional quantum Hall systems on a 
disk at filling fraction v = 5/2, with Coulomb interaction and background confining potential. We 
find the Moore-Read ground state is stable within a finite but narrow window in parameter space. 
The corresponding low-energy excitations contain a fermionic branch and a bosonic branch, with 
widely different velocities. A short-range repulsive potential can stabilize a charge +e/4 quasihole at 
the center, leading to a different edge excitation spectrum due to the change of boundary conditions 
for Majorana fermions, clearly indicating the non-Abelian nature of the quasihole. 
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Fractional quantum Hall (FQH) liquids represent 
novel states of matter with non-trivial topological or- 
der whose consequences include chiral edge excita- 
tions and fractionally charged bulk quasiparticles that 
obey Abelian or non-Abelian fractional statistics. It has 
been proposed that the non-Abelian quasiparticles can be 
used for quantum information storage and processing in 
an intrinsically fault-tolerant fashion 0, Q , in which in- 
formation is stored by the degenerate ground states in the 
presence of these non-Abelian quasiparticles, and unitary 
transformations in this Hilbert space can be performed 
by braiding the quasiparticles 0, HJ . While many Abelian 
FQH states have been observed and studied in detail [lj, 
thus far there have been relatively few candidates for the 
non-Abelian ones. The most promising candidate is the 
FQH state at Landau level filling fraction u = 5/2@. 
The leading candidate for the ground state of this sys- 
tem is the Moore- Read (MR) paired state 0], which has 
been shown 0, Q to support fractionally charged, non- 
Abelian quasiparticles. The MR state received strong 
support from numerical studies using sphere or torus ge- 
ometries @. It has been proposed that the non-Abelian 
nature of the quasiparticles in the MR state may be de- 
tected through interference experiments in edge trans- 
port 1^, U | • In order to have a quantitative understand- 
ing of the edge physics however, one needs to study the 
interplay between electron-electron interaction and con- 
fining potential, which may lead to edge structures that 
are more complicated than those predicted by the sim- 



plest theory [12J. This was found to be the case rather 



13] 



generically in the FQH regime 

In anticipation of experimental studies, in this work we 
perform detailed numerical studies of edge excitations in 
the 5/2 FQH state in finite-size systems with disc geome- 
try, taking into account the inter-electron Coulomb inter- 
action and a semi-realistic model of the confining poten- 
tial due to neutralizing background charge. For a limited 
parameter space, we find the ground state has substantial 
overlap with the MR state. Within this parameter space 
we identify the existence of chiral fermionic and bosonic 



edge modes, in agreement with previous prediction. We 
find the fermionic mode velocity is much lower than that 
of the bosonic mode. With suitable short-range repul- 
sive potential at the center, we show that a charge +e/4 
quasihole can be localized at the center of the system, 
and its presence changes the spectrum of the fermionic 
edge mode. This confirms the existence and non-Abelian 
nature of such fractionally charged quasiparticles. 

The microscopic model. We consider a microscopic 
model of a two-dimensional electron gas (2DEG) confined 
to a two-dimensional disk, with neutralizing background 
charge distributed uniformly on a parallel disk of radius 
a at a distance d above the 2DEG. This distance param- 
eterizes the strength of the confining potential, which 
decreases with increasing d. For v — 5/2, we explicitly 
keep the electronic states in the first Landau level (ILL) 
only, while neglecting the spin up and down electrons in 
the lowest Landau level (OLL), assuming they are inert. 
The amount of positive background charge is chosen to 
be equal to that of the half-filled ILL, so the system is 
neutral. The choice of a = y/AN, in units of Ib (mag- 
netic length), guarantees that the disk encloses exactly 
2N magnetic flux quanta for N = 2P, corresponding 
to v = 1/2 in the ILL 14]. The rotationally invariant 
confining potential comes from the Coulomb attraction 
between the background charge and the electrons. Using 
the symmetric gauge, we can write down the following 
Hamiltonian for the electrons confined to the ILL: 



H C = \ V mn c L+l4t c n+lCm + 



E 



U m c, n , c r . 



(1) 



where cj„ is the electron creation operator for the ILL 
single electron state with angular momentum m, V^'s 
are the corresponding matrix elements of Coulomb inter- 
action for the symmetric gauge, and C/ m 's are the matrix 
elements of the confining potential. 

MR ground state. We diagonalize the Hamiltonian 
[Eq. {!])] for each Hilbert subspace with total angular mo- 
mentum M, and obtain the ground state energy E(M). 
Figure HJa) shows E(M) vs M for N = 12 electrons 
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FIG. 1: (Color online) (a) Ground state energy E(M) in 
each angular momentum M subspace for N = 12 electrons 
with 22 orbitals in the ILL (corresponding to v — 1/2). The 
positive background charge is at a distance d = 0.4, 0.6, and 
0.8, in units of Zb, above the electron plane. The total an- 
gular momentum of the global ground state M m (indicated 
by arrows) increases from 121, 126, to 136, respectively. The 
global ground state at d = 0.6 has the same total angular 
momentum Mmr = N(2N — 3)/2 = 126 as the corresponding 
TV = 12 MR state. The overlap between the two states is 0.47. 
The curves for d = 0.4 and 0.8 have been shifted verically by 
0.6 and -0.6, respectively, (b) The total angular momentum 
of the global ground state M gs as a function of d. The arrow 
indicates the plateau at which M gs = Mmr. 



in the ILL in 22 orbitals, which is the minimum to ac- 
commodate the corresponding MR state. As illustrated 
in Fig. 01b), M gs increases with increasin g d . This is 
very similar to what happens in the 0LL [13|, and re- 
flects the interplay between electron-electron Coulomb 
repulsion and confining potential; as the confining po- 
tential weakens with increasing d, electrons tend to move 
outward, resulting in bigger M gs . M gs coincides with 
the total angular momentum of the N = 12 MR state 
M MR = N(2N - 3)/2 = 126 only within a small win- 
dow: 0.51 < d/l-Q < 0.76. This contrasts with the situa- 
tion for a Laughlin filling fraction v = 1/3, where M gs — 
N(N — l)/2v (same as the corresponding Laughlin state) 
for a substantially bigger window d < d c » 1.51b [H|. 
When the global ground state has the same total angu- 
lar momentum as the MR state, the overlap K^gsl^MR)! 2 
between the two is about 0.47. We note that, in the ab- 
sence of confining potential, the corresponding overlaps 
are 0.46 and 0.45 for N — 12 and 14, respectively. These 
values are quite substantial, given that the size of Hilbert 
subspaces are 16,660 and 194,668. But they are well be- 
low the overlap (> 0.95) in the case of the Laughlin filling 

1/3 m. The reduced overlap and window for the 
MR state reflect the fact that the paired state is much 
weaker compared to the Laughlin state. The overlap is, 
however, quite sensitive to small changes of system pa- 
rameters; for example, we can increase the overlap to 
above 0.7 for N = 12, by choosing a = \/4iV — 4, and 
a change in the V\ pseudopotential, SV\ — 0.03. Such 
sensitivity suggests that the MR state is rather "fragile" , 
consistent with experiments at v = 5/2. 



Edge excitations. The MR state has non-trivial topo- 
logical order. While our numerical results indicate that 
the ground state has substantial overlap with the MR 
state for properly chosen system parameters, it does not 
directly reflect the topological order of the system. One 
way to probe the topological order is to study edge ex- 
citations, which is also of vital experimental importance. 
For comparison, the Laughlin state supports one bosonic 
branch of chiral edge excitations, whose properties have 
been studied in tunneling experiments [15j. For v = 5/2, 
a neutral fermionic branch of excitations has been pre- 
dicted in addition to a bosonic branch [l], [l6[ . The exis- 
tence of both branches makes the low-energy excitation 
spectrum of a microscopic model at v — 5/2 richer, and 
their experimental consequences more interesting 171 ] . 

Figure [U[a) shows the low-energy excitations for pure 
Coulomb interaction and the confining potential with 
d = 0.6 for 12 electrons in 26 orbitals. Apparently, there 
is no clear distinction between fermionic and bosonic edge 
modes as well as bulk modes, due to the relatively small 
bulk gap, and system size. The situation here is similar 
to a related study on a rotating Bose gas [Hj , and will be 
analyzed in detail elsewhere. Here we focus instead on a 
model with mixed Coulomb interaction and 3-body inter- 
action for clarity. The 3-body interaction alone generates 
the MR state as its exact ground state with the smallest 
total angular momentum. The mixed Hamiltonian is 
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(1 - X)H C + XH; 
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[VfV^-r,-)^ 



(2) 

r*)], (3) 
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where S is a symmetrizer: S 123 [f 123} = A23 + /231 + 
/3i2- We measure energies in units of e 2 /els- As we 
will see, the mixed interaction, which enhances the bulk 
excitation gap with respect to that of edge excitations, 
allows for a clear separation between the two in a finite 
system, effectively increasing the system size. 

Figure HJb) shows the low-energy excitations 
AE(AM) for 12 electrons in 26 orbitals in the ILL for 
the mixed Hamiltonian with A = 0.5 and d = 0.6/b- 
There is a clear separation of the spectrum around 
AE = 0.1, below which we identify as edge modes. The 
total numbers of these states are 1, 1, 3, 5, and 10 for 
AM = 0-4, which agree with the numbers of edge states 
expected for the MR state [III]. Notably, the lowest two 
levels for AM — 4 lie very close to each other. 

In this case, we can further separate the fermionic 
and bosonic branches of the edge states. The proce- 
dure is similar to but more complicated than the one 
we used [l3| to identify edge modes in the Laughlin case 
at v — 1/3, where there is only one bosonic branch of 
edge modes. The basic idea is to label the low-lying 
states by two sets of occupation numbers {nb(Zb)} and 
{rif(lf)} for bosonic and fermionic modes with angular 
momentum ?b,fj respectively. Since the fermionic edge 
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FIG. 2: (Color online) Low-energy excitations AE(AM) 
from exact diagonalization (solid lines) for N = 12 electrons 
in 26 orbitals in the ILL (corresponding to v = 1/2) for (a) 
Coulomb Hamiltonian [Eq. JT}] and (b) mixed Hamiltonian 
[Eq. ((2| with A = 0.5]. The neutralizing background charge 
for the Coulomb part is deposited at d — 0.61b above the 
electron plane, (c) Dispersion curves of bosonic (J5b) and 
fermionic modes (Ef) of the system. These energies can be 
used to construct the complete edge spectrum for the 12- 
electron system up to AM = 4 [dashed bars in (b)]. 



excitations are Majorana fermions that obey antiperi- 
odic boundary conditions [lq |. k must be positive half 
integers and ^\ n f(k) an even integer because fermion 
modes are occupied in pairs. The angular momentum 
and energy of the state, measured from those of the 
ground state, are AM = J2i h n b(lb)lb + J2i { n f('f)'f and 
A-E = Zi b n h (l h )E h (l h ) + J2 h n t (k)E { (k). The difficulty 
here is, in addition to the bosonic modes, we also have 
fermionic modes, and thus the convolution of fermionic 
and bosonic modes. Fortunately, we note that the num- 
ber of states with nearly zero energy coincides with the 
number of fermionic edge states expected by theory. Ac- 
cordingly in our construction we will assume these ener- 
gies to have been evolved from combining two Majorana 
fermions. Through careful analysis of the low-energy ex- 
citations up to AM = 4, we obtain the results of bosonic 
and fermionic mode energies up to lb = 4 and If = 7/2, 
respectively, plotted in Fig. HJc). The detailed analysis 
will be published elsewhere. Using these 8 energies [ex- 
cluding the trivial £b(0) = 0], we can construct the whole 
low-energy spectrum of the system up to AM = 4, a to- 
tal of 20 states. The excellent agreement [see Fig. HJb)] 
justifies our analysis, and thus our result, that, energet- 
ically, fermionic modes are well separated from bosonic 
modes. In contrast to the roughly linear dispersion of the 
fermionic branch, the energy of the bosonic branch bends 
down (despite a much bigger initial slope or higher veloc- 
ity), suggesting a potential vulnerability to edge recon- 
struction in the bosonic branch 



13j. These are not sur- 



prising since the bosonic modes are charged; as a result 
its velocity is dominated by the long-range nature of the 



Coulomb interaction in the long- wavelength limit, but in 
the meantime it is also more sensitive to the competi- 
tion between Coulomb interaction and confining poten- 
tial which can lead to instability at shorter wavelength. 

Charge +e/4 and +e/2 quasiholes. One of the most 
important properties of the MR state is that it sup- 
ports charge ±e/4 quasihole/particle excitations. To 
demonstrate the unusual fractional charge, we add, to 
the mixed Hamiltonian with A = 0.5 and d = 0.5/b, 
a short-range potential: H\y — Wc\cq, which tends to 
create quasiparticles or quasiholes at the origin. For 
small enough repulsive W, the ground state of the sys- 
tem should remain MR-like. As W is increased, a quasi- 
hole of charge +e/4 can appear at the origin, reflected 
by a change of ground state angular momentum from 
M gs = N(2N - 3)/2 to N(2N - 3)/2 + N/2, and de- 
pletion of 1/4 in the total occupation number of elec- 
trons at orbitals with small angular momenta. If W is 
increased further, a +e/2 quasihole, much like a quasi- 
hole for the Laughlin state, appears near the origin in 
the global ground state, whose total angular momentum 
further increases to N (2N — 3) /2 + N . This is observed 
for a system of 12 electrons in 24 orbitals (as well as a 
smaller system of 10 electrons in 20 orbitals). Figure[3ja) 
shows the increase of M gs from 126 to 132 and then to 
138 with increasing W. Fig. [31(b) compares the electron 
occupation number n(m) in each orbital for W = 0.0 
(M^ = 126, MR-like) and W = 0.1 (M gs = 132). The 
accumulated difference in the occupation numbers of the 
two states, X^o^ n (*)' oscillates around —0.25 for rn 
up to about 19, indicating the existence of a +e/4 quasi- 
hole at the origin. The same comparison for W = 0.1 
(M gs = 132) and W = 0.25 (M gs = 138) is plotted in 
Fig. [3jc). Their difference (~ —0.25) indicates the emer- 
gence of another +e/4 quasihole at the origin, or a +e/2 
quasihole compared to the MR-like state for W = 0.0. 

The +e/4 quasihole supports a zero energy Majorana 
fermion mode which is responsible for its non-Abelian na- 
ture. This zero mode pairs with the edge excitations and 
changes their spectra. With quasiholes in the bulk, the 
fermionic edge excitations are Majorana fermions that 
obey either periodic (integer If, twisted sector) or an- 
tiperiodic (half integer If, untwisted sector) boundary 
conditions (l6| . In the presence of an odd number of 
+e/4 quasiholes (twisted sector), two- fermion edge exci- 
tations are shifted by S (AM) = — 1, relative to those in 
the presence of an even number of such quasiholes (un- 
twisted sector). This is demonstrated in Fig. fS)Jd)-(f) for 
12 electrons in 24 orbitals. From A/ gs = 126 in Fig. EJd) 
(W = 0.0, no quasihole present), we count the numbers of 
fermionic edge states as 0, 1, 1, 2 for AM — 1-4. From 
M gs = 132 in Fig. [Sfe) (W = 0.1, one +e/4 quasihole 
present), the numbers change to 1, 1, 2, 2 for AM = 1- 
4. In particular, the existence of a fermionic state at 
AM = 1 clearly indicates the change. From M gs = 138 
in Fig. [3{f) (W — 0.25, two +e/4 quasiholes present), 
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FIG. 3: (Color online) Generation of quasiholes using a 
short-range repulsion VFcJco in a system of 12 electrons in 
24 orbitals, for the mixed Hamiltonian Eq. ([2]) with A = 0.5 
and d = 0.5/b. (a) Ground state angular momentum M gs 
as a function of W . (b) Electron occupation number n(m) 
of the ground state for W = 0.0 (M gs = 126) and W = 0.1 
(M ga = 132), as well as the accumulated difference in n(m) 
between the two states, X^o ^ n (*)i which oscillates around 
-0.25 (dotted line) for m up to about 19, indicating the emer- 
gence of a charge +e/4 quasihole at m = 0. (c) n(m) for 
W = 0.1 and W = 0.25 (M gs = 138), and their accumu- 
lated difference, indicating the appearance of another +e/4 
quasihole. Low-energy edge excitations of the systems are 
plotted for (d) W = 0.0, (e) W = 0.1, and (f) W = 0.25. 
The fermionic mode supports 0, 1, 1, 2 states (solid bars) for 
AM = 1-4 in (d) (M gs = 126), but 1, 1, 2, 2 for AM = 1-4 
in (e) (M gs = 132). The numbers change back to 0, 1, 1, 2 
for AM = 1-4 in (f) (M gs = 138). This suggests that a sin- 
gle +e/4 quasihole (or an odd number of quasiholes) changes 
the fermionic mode spectrum while a +e/2 quasihole (or, in 
general, an even number of +e/4 quasiholes) does not. 



the numbers change back to 0, 1, 1, 2 for AM = 1-4. 

Summary. Our results suggest that the Moore-Read 
(MR) state properly describes a half-filled first Landau 
level, for properly chosen confinement potential. In this 
case the system supports chiral edge excitations as well 
as fractionally charged quaisholes, and their properties 



agree with theory predictions. We also find that the win- 
dow of stability of the MR state is rather narrow, and 
the edge modes may suffer from reconstruction or other 
instabilities as the confinement potential varies. The na- 
ture and consequences of such instabilities are currently 
under investigation, which will be presented elsewhere 
along with further details of the present work. 
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